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Abstract 

The Dirac equation with a scalar and an electromagnetic potentials is considered. In the time- 
harmonic case and when all the involved functions depend only on two spatial variables it reduces 
to a pair of decoupled bicomplex Vekua-type equations [§]. Using the technique developed for 
complex Vekua equations a system of exact solutions for the bicomplex equation is conctructed 
under additional conditions, in particular when the electromagnetic potential is absent and the 
scalar potential is a function of one Cartesian variable. Introducing a transmutation operator 
relating the involved bicomplex Vekua equation with the Cauchy-Riemann equation we prove 
the expansion and the Runge approximation theorems corresponding to the constructed family 
of solutions. 
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1 Introduction 

The Dirac system with a scalar and an electromagnetic potentials is considered. In [8] (see also [18] ) 
it was shown that in the time-harmonic case and when the whole model is independent of one of 
the spatial variables the system reduces to a pair of decoupled Vekua-type equations which differ 
from the classical Vekua equations considered in the theory of generalized analytic or pseudoanalytic 
functions [2J, [TS], [3S] by the fact that they are bicomplex. In 8 using this reduction as well as a 
procedure introduced by L. Bers, for an arbitrary scalar potential depending on one Cartesian variable 
an infinite family of solutions of the Dirac system was constructed. Nevertheless the completeness of 
this family in the linear space of all solutions was not proved due to the lack of some fundamental 
results in the theory of bicomplex Vekua equations such as the similarity principle and many other. 

The constructed family of solutions is a system of formal powers generalizing those introduced 
by L. Bers onto the bicomplex situation. Meanwhile in the classical complex case there is a well 
developed theory of formal powers with the Runge-type approximation theorem and other related 
results (see [5] and references therein), in the bicomplex case up to now no such result was available 
even for simplest examples. 

In the present work we prove the completeness of the family of solutions obtained in [5] by using 
so-called transmutation operators and some recent results on their mapping properties [B.. The 
notion and the name of the transmutation operator appeared in the work of J. Delsarte [llj . [12] and 
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later developed in [T3J, [23J and many other publications [T], [7J, [22], [23], [25], [27J. Combining the 
results from [6] on mapping properties of the transmutation operators with the results from [21] on the 
construction of a transmutation operator for the Darboux transformed Schrodinger equation we obtain 
transmutation operators which relate the bicomplex Vekua equations arising from the Dirac system 
with the Cauchy-Riemann equation. Using this result we prove that the bicomplex pseudoanalytic 
formal powers are the result of application of a corresponding transmutation operator to the usual 
powers of the complex variable z. This together with the boundedness of the transmutation operator 
and of its inverse allows us to prove the expansion and the Runge approximation theorems for solutions 
of the considered bicomplex Vekua equations. 



2 The Dirac system and bicomplex pseudoanalytic functions 

Consider the Dirac operator with a scalar and an electromagnetic potentials 

E> = 7o 9 t + X! 1k d k + i I m + p e a + ^2 A klk + Psc j 
fc=i V fe=i / 

where 7-, j = 0, 1, 2, 3 are usual 7-matrices (see, e.g., [4J, [26] ) 
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to G M, pel, Ak and p sc are real valued functions. 

We will denote the algebra of biquaternions or complex quaternions by H(C) with the standard 
basic quaternionic units denoted by eo = 1, e±, e% and e%. The complex imaginary unit is denoted by 
i as usual. The set of purely vectorial quaternions q — q is identified with the set of three-dimensional 
vectors. 

The quaternionic conjugation of a biquaternion q = qo + q will be denoted as q = qo — q. 
Sometimes the following notation for the operator of multiplication from the right-hand side will be 
used M p q — q ■ p. 

The main quaternionic differential operator introduced by Hamilton himself and sometimes called 
the Moisil-Theodoresco operator is defined on continuously differentiable biquaternion-valued func- 
tions of the real variables X\, x% and x 3 according to the rule 

3 

Dq = ^2 e kdkq, 

k=l 



where d, 




In [14] (see also [9], [15], [20]) a simple invertible matrix transformation was obtained which allows 
one to rewrite the classical Dirac equation in biquaternionic terms. Namely, the Dirac operator D is 
equivalent to the biquaternionic operator 

R = D — d t M ei + a + M- i < ? " ei - i ®» +m > e3) 
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where a = i(A%ei + A2e 2 — ^.363) and the notation means the reflection with respect to X3, 

f := f(t, x\, X2, — X3). Note that in the absence of the electromagnetic potential the operator R 
becomes real quaternionic which is an important property (see |19j). 

In what follows we assume that potentials are time-independent and consider solutions with a 
fixed energy: <!>(£, x) = $ w (x)e lwt . The equation for $ w has the form 



®A = in G (1) 

where G is a domain in R 3 , 

D w = iuj-/ a + ^2 7aA- + i I m + p e r/ + A k J k + p sc ) ■ 

fc=l \ k=l I 

Under the mentioned above matrix transformation the operator B> u turns into its biquaternionic 
counterpart 

R u = D + a + M b 

with b = —i((p e i + uS)ei — i(p sc + rnje-z). Thus, equation ([T|) turns into the complex quaternionic 
equation 

R^q = (2) 

where q is a complex quaternion valued function. In what follows we study this equation. 

Let us introduce the following notation. For any biquaternion q we denote by Qi and Q2 its 
bicomplex components: 

Qi=qa+q3e 3 and Q 2 = <?2 - <7ie 3 - 

Then q can be represented as follows q = Qi + Q2&2- For the operator D we have D = D\ + D^e-i 
with D\ — e-idz and D2 = 82 — <9ie3. Notice that b = Be2 with B = —(p sc + m) + i(p e i + w)e3, 
a = A\ + A2e2 with A\ = 0363 and A2 = 0,2 — axes- 

We obtain that equation ((2J is equivalent to the system 

D 1 Q 1 -D 2 Q2 + A 1 Q 1 -A2Q2-BQ2 = 0, (3) 

D 2 Q X + D1Q2 + A 2 Q X + A1Q2 + BQx = 0, (4) 

where Q\ and Q2 are bicomplex components of q. We stress that the system ((3j) , ((4j) is equivalent 
to the Dirac equation in 7-matrices (fT|). 

Let us suppose all fields in our model to be independent of X3, and A\ = 0.363 = 0. Then the 
system ^ decouples, and we obtain two separate bicomplex equations [5], [TB] 

D2Q2 = -A2Q2 - BQ 2 , and ^ 2 Qi = -A2Q1 ~ ~BQx- 

Denote d = D2, a = —A2, b = —B, w — Q2, W = Q±, z = x + yk, where x — X2, y = X\ and for 
convenience we denote k = e^. Then we reduce the Dirac equation with electromagnetic and scalar 
potentials independent of X3 to a pair of Vekua-type equations 

dw = aw + bw (5) 

and 

dW = aW + bW. (6) 
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3 Some definitions and results from bicomplex pseudoana- 
lytic function theory 

Definition 1 We consider M -valued functions of two real variables x andy. Denote d — +k^-) 



dy i 

dy> 



and d = \{-§^ — kj^-). An equation of the ft 



dw = aw + bw, (7) 

where w, a and b are M-valued functions is called a bicomplex Vekua equation. When all the involved 
functions have their values in Ck only, equation £?J) becomes the well known complex Vekua equation 
(see U8f . J28f ). We will assume that w G C (fi) where SI C I 2 is an open domain and a, b are 
Holder continuous in f2. 

When a = and b = ^ where 4> : — > Q possesses Holder continuous partial derivatives in il 
and <f>{x, y) ^ 0, V(x, y) £ £1 we will say that the bicomplex Vekua equation 

dw — ^-w (8) 

<P 

is a Vekua equation of the main type or the main Vekua equation. 

For classical complex Vekua equations Bers introduced [2] the notions of a generating pair, gen- 
erating sequence, formal powers and Taylor series in formal powers. As was shown in [5], |18j the 
definition of these notions can be extended onto the bicomplex situation. Here we briefly recall the 
main definitions. 

Definition 2 A pair of M-valued functions F and G possessing Holder continuous partial derivatives 
in fl with respect to the real variables x and y is said to be a generating pair if it satisfies the inequality 

Vec(FG) ^ in Q. (9) 

Condition (|9]) implies that every bicomplex function w defined in a subdomain of f2 admits the 
unique representation w = <pF + %pG where the functions <fi and are scalar (C^-valued) . 

Remark 3 When F = 1 and G = k the corresponding bicomplex Vekua equation is 

dw = 0, (10) 

and its study in fact reduces to the complex analytic function theory Indeed, consider the following 
pair of idempotents P + = |(1 + ik) and P = ^(1 — ik) ((P^) = P^)- Then the functions P + w 
and P~w are necessarily antiholomorphic and holomorphic respectively. Indeed, application of P + 
and P to \10\) gives us 

d z P + w = and i\P^w = (11) 

where d z = ^(^~i~§y~) andi\ = \(~§^+i-~§^)- Moreover, P + w = P + (u+jv) = P + (u~ iv) andP~w = 
P (u + iv) . Due to Ml)) the scalar functions w + := u — iv and w^ := u + iv are antiholomorphic and 
holomorphic respectively. We stress that w + is not necessarily a complex conjugate of w~ (u and v 
are Ci-valued). 

Let us notice that due to the equivalence of \10\) and ill]) we have that a bicomplex solution w of 
\10\) admits a convergent Taylor series w(z) = Y^=o a « z ™ */ an< ^ only if the series X^^Lo a n ( z+ ) 
and X)^Lo a n ( z ~) n corresponding to w + and w~ respectively converge (here and are scalars, 
a„ = Sc(a n ) =F iVec(a n ) and z^ = x =p iy). In particular, the radius of convergence of the power 
series X)^Lo a ™ z ™ ^ as ^ e f orm R = min{i? + ,i?_} where 1/R± = limsup JWOO la^j 1 ^™. 

Assume that (F, G) is a generating pair in a domain f2. 
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Definition 4 Let the M -valued function w be defined in a neighborhood of zq G fl C Ck- In a complete 
analogy with the complex case we say that at zq the function w possesses the (F,G) -derivative w{zq) 
if the (finite) limit 

w(zo) = Hm *W ~ X ° F ^ - (12) 

Z^rZo Z — Zq 

exists where Ao and /i are the unique scalar constants such that w(zq) = XqF(zq) + [i G(zo). 

Similarly to the complex case (see, e.g., jTSJ Chapter 2]) it is easy to show that if w(zq) exists 
then at Zq, dw and dw exist and equations 

dw = a( Ft G)W + b(F t G)W (13) 

and 

w = dw — Afp^Q-jW — Brp^ayw (14) 

hold, where cL(f,g)j ^(f.g)j ^(f,g) an d B(f,g) are the characteristic coefficients of the pair (F,G) 
defined by the formulas 

FdG-GdF , FdG-GdF 



a{F,G) == — ==t^, b{F,G) 



FG-FG y ' ' FG-FG 

. FdG-GdF FdG-GdF 

A(F > G) ~ ~ FG-FG ' B(F > G) ~ FG-FG ' 

Note that FG-FG= -2k Vec(FG) ^ 0. 

If dw and dw exist and are continuous in some neighborhood of Zq, and if (|13j) holds at zq, then 

w{zq) exists, and (fT4| holds. Let us notice that F and G possess (F, G)-derivatives, F = G = and 
the following equalities are valid which determine the characteristic coefficients uniquely 

dF = a(F,G)F + b(F,G)F, dG = a (F . G) G + b( F . G) G, 

OF = A {F>G) F + B(p G yF , dG = A (FjG) G + B {F , G) G. 

If the (F, G)-derivative of a B- valued function w — <fiF + ipG (where the functions <f> and ip are scalar) 
exists, besides the form flT4l it can also be written as follows w = dcj) F + dip G. 

Definition 5 Let (F,G) and (iq,Gi) - be two generating pairs in tt. (iq,Gi) is called successor 
of (F, G) and {F, G) is called predecessor of (iq, Gi) if 

a (F 1 ,G 1 ) = a (F,G) and &(Fi,Gi) = ~B(F,G)- 
By analogy with the complex case we have the following statement. 

Theorem 6 Let w be a bicomplex {F,G)-pseudoanalytic function and let (iq,Gi) be a successor of 
(F,G). Then w is a bicomplex (Fi,G\)-pseudoanalytic function. 

Definition 7 Let (iq G) be a generating pair. Its adjoint generating pair (F,G)* = (F*,G*) is 
defined by the formulas 

2i? 2G 

F = = — G = — = — 

FG FG FG-FG 
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The (F, G)-integral is defined as follows 

J Wd {FtG) z=^F(z 1 )Sc J G*Wdz + G{ Zl )Sc J F*Wdz^j 

where T is a rectifiable curve leading from zq to z\. 

If W — <pF + ipG is a bicomplcx (F, G)-pseudoanalytic function where 4> and ip are complex valued 
functions then 

Wd (F , G) z = W(z) - <P{z )F{z) - i>{z )G{z), (15) 



/ 

•> Zn 



and this integral is path-independent and represents the (F, G)-antiderivative of W. 

Definition 8 A sequence of generating pairs {(F m , G m )}, in = 0, ±1, ±2, . . . , is called a generating 
sequence if (F m+1 ,G m+ i) is a successor of (F m ,G m ). If (Fo,Go) = (F,G), we say that (F,G) is 
embedded in {(F m ,G m )}. 

Let If be a bicomplex (F, G)-pseudoanalytic function. Using a generating sequence in which 
(F, G) is embedded we can define the higher derivatives of W by the recursion formula 

W M=W; W ^\ = d (^ 7 ) W[m \ m= i, 2 ,.... 

dz 

Definition 9 The formal power Zm\a, Zq\ z) with center at zq G £1, coefficient a and exponent 
is defined as the linear combination of the generators F m , G m with scalar constant coefficients X, jj, 
chosen so that XF m (z ) + fiG m (z ) — a. The formal powers with exponents n = 0, 1, 2, . . . are defined 
by the recursion formula 

Z( n " +1 >(a,z ;z) = (n+1) f Z^ +1 (a,z ;Od {Fm ,G m )C (16) 

J z 

This definition implies the following properties. 

1. Zm\a, z ; z) is an (F m , G m )-pseudoanalytic function of z. 

2. If a' and a" are scalar constants, then 

ZW (a' + ka", z 0] z) = a'zW (1, z ; z) + a"Z™ (k, z ;z). 

3. The formal powers satisfy the differential relations 

d(F m ,G m )Z ( rn\a,Z ;z) _ „(„-!), . > 
-J~ z ~~ nZ, rn+l 1°' Z 0, Z)- 

4. The asymptotic formulas 

Z^\a,z ;z) ~ a(z - z ) n , z -> z 

hold. 

The case of the main bicomplex Vekua equation is of a special interest also due to the following 
relation with the stationary Schrodinger equation. 
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Theorem 10 U6]/ Let W = W\ + lsW 2 be a solution of the main bicomplex Vekua equation 



dW = in CI (17) 

where W\ — ScW , W 2 = VecW and the <Ci-valued function <j) is a nonvanishing solution of the 
equation 

— Au + qi(x, y)u = in CI (18) 

where q\ is a continuous Ci-valued function. Then W% is a solution of U8\) in Cl and W 2 is a solution 
of the associated Schrddinger equation 

- Av + q 2 (x 7 y)v = in Cl (19) 

where q 2 = 8^^- - qi . 

We need the following notation. Let w be a B-valued function defined on a simply connected 
domain Q with W\ — Sew and w 2 = Vec w such that 

and let T C fl be a rectifiable curve leading from (xo,yo) to (x, y). Then the integral 

Aw(x, y) := 2 (^J w\dx + w 2 dy^j 

is path- independent, and all C^-valued solutions ip of the equation dtp = w in 57 have the form 
ip(x,y) — Aw(x,y) + c where c is an arbitrary C^-constant. In other words the operator A denotes 
the well known operation for reconstructing the potential function from its gradient. 

Theorem 11 flffj Let W\ be a Ci-valued solution of the Schrddinger equation A18\) in a simply 
connected domain fl. Then a Ci-valued solution W 2 of the associated Schrddinger equation U9\) such 
that Wi + \lW 2 is a solution of |-?7| ) in CI can be constructed according to the formula 

W,= l -A(^d(^-\\ + ^ 



where c\ is an arbitrary Ci-constant. 

Vice versa, given a solution W 2 of H9\) , the corresponding solution W% of H8\) such that Wj+kW-j 
is a solution of |17[ ) has the form 

Wx = -<j>A(^d(4>W 2 ) S J +c 2 

where c 2 is an arbitrary Ci-constant. 

As was shown in [17] (see also [18]) a generating sequence can be obtained in a closed form, for 
example, in the case when <p has a separable form (f> = S(s)T(t) where s and t are conjugate harmonic 
functions and S, T are arbitrary twice continuously diffcrcntiablc functions. In practical terms this 
means that whenever the Schrddinger equation (1181) admits a particular nonvanishing solution having 
the form (f> — /(£) g(rj) where (£, 77) is one of the encountered in physics orthogonal coordinate systems 
in the plane a generating sequence corresponding to (fl7j) can be obtained explicitly [18] Sect. 4.8] . The 
knowledge of a generating sequence allows one to construct the formal powers following Definition [9] 
This construction is a simple algorithm which can be quite easily and efficiently realized numerically 
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[5], [TO]. Moreover, in the case of a complex main Vekua equation which in the notations admitted in 
the present paper corresponds to the case of <f> being a real- valued function (then the main bicomplex 
Vekua equation decouples into two main complex Vekua equations) the completeness of the system 
of formal powers was proved [5] in the sense that any pseudoanalytic in fi and Holder continuous 
on <9£! function can be approximated uniformly and arbitrarily closely by a finite linear combination 
of the formal powers. The real parts of the complex pseudoanalytic formal powers represent then a 
complete system of solutions of one Schrodinger equation meanwhile the imaginary parts give us a 
complete system of solutions of the associated Schrodinger equation. 

4 Transmutation operators and a complete family of solutions 
of the Dirac equation 

Let us consider the following situation p sc — p(x), p e i — and ^4=0. Then the Dirac equation is 
equivalent to the pair of bicomplex Vekua equations 

dw = bw (21) 

dW = W (22) 

where b = p(x) + m — zwk. 

Let P denote an antiderivative of p. Consider the function 

4>(x, y) = e p (*)+™+^y 

Note that dejy/cf) = b. Then if W is a solution of (f2"2l then the complex valued function W\ = Sc(W) 
is a solution of the Schrodinger equation 

(-A + v) Wi = 0, with v(x) =p'(x) + (p(x) + m) 2 - J 2 

and the complex valued function Wi = Vec(W) is a solution of the associated Schrodinger equation 

(-A + fi) W 2 = 0, with (i{x) = -p'(x) + (p(x) + mf - uj 2 

On the other hand equation (|22l) can be written as a main Vekua equation 

9-^cV = (23) 



where 

4>{x,y) = f(x)g(y) with f(x)=e p ^ +mx and g(y) = 

Notice that / and g are complex valued functions. We assume that their domains of definitions are 
finite segments [—a, a] and [—6, b] respectively. Assuming that p £ C 1 [—a, a] we obtain that / and 
g are nonvanishing C 2 -functions. The separable form of <j> allows us to write down a generating pair 
associated with equation (|23|) (F.G) — (</>, k/0) as well as a generating sequence of the period two 
embedding this generating pair 

(F,G) = (<t>,k/ct>); (F u G 1 ) = (0/f 2 ,kf 2 /4>) 
(F 2 ,G 2 ) = (F, G) ■ (F 3 ,G 3 ) = (F 1 ,G 1 ) ; .... 

The corresponding formal powers can be constructed as follows. We consider the formal powers with 
the centre in the origin and for simplicity assume that /(0) = 1 (for g this is also the case). Define 
the following systems of functions {<fik}T=o an< ^ Wfclfclo 

, . r f{x)x( k \ x ), k odd 

lf{x)X\ K ){x), k even 



where 



= x (0 Hx) = i 

X^(x) = nf ^(p)[f»W] ( - ir * 
xM(x) = n T X^){p)[f{pf- ir ~ X dp 



and 



where 



( ((/|= .g(y)YW(y),k odd 
\g(y)Y( k \y), k even 



y(°)( y ) = r(°)( y ) = i 

rW(y) = nf y<»-%)[^(0] ( - 1)B_I de. 



Then the formal powers based on the given generating sequence are defined by the formulas 

(a, 0; z) = c(>[x, y)Sc*Z^ (a, 0; z) + ,, k ^ec.Z^ (a, 0; z) (26) 



where 



a 'H (™)x("- m )k m y( m ) +ka'^] (»)x("- m )k ro y( m ), n odd 
,zW(a,0;«) = | ™=° ™=° ( 27 ) 

(^)X("- m )k m y( m ) +ka ^ (^)X("- m )k m y( m ), n even. 

m— m— 

In a similar way the formal powers corresponding to (|2ip can be constructed, they will be denoted 
as Z (a, 0; z). Notice that a generating pair for (|2"Tj) is given by 

~ 9 ~ / 

F = - and G = k-. 

/ 5 

In [B] it was shown that for functions / and g satisfying the above conditions there exist the 
transmutation operators Tf and T g defined as follows 



T f [u(x)\ = u(x) + / K{x,t;f(0)).u(t)dt (28) 



where the kernel K(x, i; /'(0)) is given by 

K(x, t; f(0)) = ^ + f) + M |* s ) - fffo -a)] 
and the function K{x,t) is the unique solution of the Goursat problem (see [24] ) 

K(x,t) = ^K(x,t) 
^K(x, x) = h Jq q(s)ds; K(x, —x) = 



and 

rV 



T g [v(y)} = v(y) + / K(y, t; g'(0))v(t)dt 

J-y 

where 

K(y, t; g(0)) = ^ + K(y, t) + ^-f [k(y, s) - K(y, -a)] ds 

and the function K(x, t) is the unique solution of the Goursat problem 

(j^-q( y ))k{y,t) = ^K{y,t) 
k(y,y) = ±f*q(s)ds; K(y,-y)=0 

with q = f"/f and q = g"/g. 

Moreover, Tf and T g satisfy the relations 

T } [x k ] = Vk and T g [y k ] = Vfc g f% (29) 

We will need similar systems of functions {^j-^Lq and < ip k > corresponding to 1// and 1/g 
respectively, 

f-J-^X^Va:), fc even 

ly^yXW^), fcodd 

(7T7iY {k) (y), k even 
^(2/)= ^- (31) 

For these systems of functions another pair of transmutations T\if and T\/ g is constructed (see 
[2 1J ) . one of the representations of which can be given by the equalities 

T 1/f u(x) = JL { f f(r})Tf [du{ri)\ d v + u(0)\ , 



m 

Ti /g v(y) = -L | jf ff fa)T a [flvfa)] d, 7 + «(0) 
They satisfy the equalities 

T v/ = £3 fc and T 1/9 [y k ] = % k , Vfc g N . (32) 
The operators Tj.// and admit the representations as Volterra integral operators |21) . 

T 1/f u(x) = u(x) + f K 2 (x,t;-f(0))u(t)dt, 

J — X 

where the kernel K2(#, £; —f(0)) has the form 

K 2 (x, t; -/t0)) = --^y ( y ° dtK^s, t; f(0))f(s) ds + ® /(-*)) 

and the formulas for 7i/ g are completely analogous with an obvious substitution of / by g. 
The introduced transmutation operators satisfy interesting commutation equalities. 



10 



Corollary 12 \21f The following operator equalities hold on C 1 -functions of the respective variables 

d x fT 1/f = fT f d x , d x jT f = jT 1/f d x . (33) 

dygTl/g = gTgdy, dy^Tg = ^ 7\ / g 8y . (34) 

Consider the operators projecting onto the scalar and the vector parts respectively 

P+ = l(I + C) and P - = ±(I-C). 

Let us introduce the following operators 

T =T f TgP+ + kT 1/f T 1/g p- (35) 

and 

Ti =Ty f T g P+ + kT f T 1/g p- . (36) 

From now on let f2 C R — [—a, a] x [—b, b) be a simply connected domain such that together with 
any point (x,y) belonging to Q the rectangle with the vertices (x,y), (—x,y), (x,—y) and (— x, — y) 
also belongs to f2. In such a domain application of operators To and Ti is meaningful. 

Proposition 13 The following equalities hold for any M-valued, continuously differentiable function 
w defined in fl. 



8 - ^-C^j T w = Ti (dw) , (d + ^-C^j T lW = T (Bw) 



(37) 



d- -jCj T w = Ti (dw), I 8 + -^-C \ Tiw = T (dw) . (38) 

Proof. The proof consists in a direct calculation with the aid of the relations from Corollary Q21 ■ 
An immediate corollary of equalities (|37j) is the fact that the operator To maps bicomplex ana- 
lytic functions into (0, k/0) — pseudoanalytic, i.e., into solutions of (1231) and the operator Ti maps 

bicomplex analytic functions into (j, -pseudoanalytic i.e., into solutions of the equation 

d+^c)w = 0. (39) 



Moreover, they map powers of the variable z into corresponding formal powers. 
Proposition 14 For any z £ fl and o£B the following equalities are valid 

T [az n ] = Z (n) (a,0;z) and Ti[az"] = z[ n) (a,0; z). 
Proof. The proof consists in the observation that for a = a'+ k&' and z = x + ky one has 



m=0 



(o'+ k60 J2 [ x n - m k m y m 



and the result follows from the formulas (f26j) . ([27| by application of the mapping properties (|29|) , 
(02]). ■ 

Notice that both To and Ti are bounded operators on the space of continuous functions with re- 
spect to the norm ||w|| = max(|u| + \v\) where w = u + kv. Indeed, consider ||Tow|| = max(|T/T s u| + 
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\TijfTi jgV\) < Mimax|w| +M2max|t;| where the constants Mi and M2 depend only on the corre- 
sponding kernels of the bounded Volterra operators Tf, T g , T x /f and Ty g . Then ||T w|| < M \\w\\ 
where M — max{Mi,M2}. The proof for the operator Ti is analogous. Moreover, Tq 1 and 1 are 
bounded as well (the form of the inverses for Tf, T g , T X jf and T\i g can be found in [3T]) due to the 
fact that their integral kernels enjoy the same regularity properties as the kernels of To and Ti. 
Let us establish another useful fact concerning the mapping properties of the operators To and 

Ti. 

Proposition 15 Let w be a bicomplex analytic function in f2 and W = Tow be a corresponding 
solution of l\23\i . Then 

T (d (2n) w) = W l2n] and T x (d^'^w) = W [2n ~ 1] , n = l,2,.... (40) 

Proof. From (f3"5| we have 

W = T x (dw) . (41) 

W is a solution of the succeeding Vekua equation (|39|). Denote W\ = W. Any solution of ([39]) is the 
image of a bicomplex analytic function under the action of the operator Ti, so W\ — T\Wi. Due to 
([38]) we have Wi = T (dwi). Thus, W = T (d 2 w) because from (gTJ) wi = dw. Now flU can be 
easily proved by induction. ■ 

The established relations from Propositions 1X31 and IT41 together with the fact that To and Ti are 
bounded operators together with their respective inverses allow us to transfer several results from 
analytic function theory onto the solutions of the bicomplex Vekua equations under consideration 
and as hence onto the solutions of the Dirac system with a scalar potential being a function of one 
Cartesian variable. Here we give two examples of such results. 

Theorem 16 Let W be a solution of (jgffp in a disk D with the center in the origin and radius R. 
Then it can be expanded into a Taylor series in formal powers 

00 

W(z) = J2z {n) (a n ,0;z) 

71=0 

with the radius of convergence R. The series converges normally in D and the coefficients a n have 
the form 

_ tyW(o) 

— ; • 

n! 

Proof. Consider w = Tq 1 W . It is a bicomplex analytic function, so we have that it can be expanded 

00 

into a Taylor series w(z) = a n z n with the coefficients a n = - /n\. Application of 'To gives us 

n=0 

00 00 

a series for W , W(z) = Tqw(z) = To[a„z n ] = (a„, 0; z). Due to the uniform boundedness 

n=0 n=0 

of To the radius of convergence of the series is preserved. Note that the Taylor coefficients coincide. In 
order to finish the proof we use Proposition \15\ and the fact that both operators To and Ti preserve the 
values of a function in the origin. This is obvious from their definition and from the Volterra integral 
form of the operators Tf, T g , T 1/f and T 1/g (see, e.g., 1%$). Thus, W^{0) = T (<9 (2n) w) (0) = 
d^wiO) and W^ n ~ l \0) = T 1 (d {2n -^w) (0) = d^-^ w{0) , u=l,2,.... ■ 

Theorem 17 Any solution W of ^23^ in VL can be approximated arbitrarily closely on any compact 

N 

subset K of il by a finite combination of formal powers (a formal polynomial) Z^{a n , 0; z). 

71=0 
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Proof. Consider w = Tg W. Due to the Runge approximation theorem the function w can be 
arbitrarily closely approximated by a polynomial in z. Then due to the boundedness of To and Tq 1 
and Proposition 1141 we obtain the required result for W . ■ 

This theorem in fact means the completeness of the family of functions 

{ZW(1,0;*), ZM(k,0;z)}°° 

I J n=0 

in the space of all solutions of the Vekua equation (|23p . A similar fact is true for equation (131)1) and 
corresponding formal powers. The combination of both families of formal powers gives us a complete 
family of solutions of the Dirac equation ([2} in the considered case. 
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